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O ■ Abstract 



We prove the analyticity of the finite volume QCD partition function for complex values of the 
^-vacuum parameter. The absence of singularities different from Lee- Yang zeros only permits A 
cusp singularities in the vacuum energy density and never V cusps. This fact together with the 
Vafa-Witten diamagnetic inequality implies the vanishing of the density of Lee- Yang zeros at i? = 
■ and has an important consequence: the absence of a first order phase transition at 6 = 0. The 

o. 

result provides a key missing link in the Vafa-Witten proof of parity symmetry conservation in 
vector-like gauge theories and follows from renormalizability, unitarity, positivity and existence of 



BPS bounds. Generalizations of this theorem to other physical systems are also discussed, with 
particular interest focused on the non-linear CP N sigma model. 



The Vafa-Witten theorem [l| is one of the very few non-perturbative analytic results 
of QCD. In absence of supersymmetry, non-Abelian continuum gauge theories in 3+1 di- 
mensions do exhibit many very interesting non-perturbative phenomena which have been so 
far elusive to analytic approaches. However, by a very simple and elegant argument Vafa 
and Witten |l( explained why we do not observe spontaneous violation of parity symmetry, 
not only in QCD but in any gauge theory with Dirac fermions in 3+1 space-time dimen- 
sions. The relevance of the physical implications of this theorem require to further clarify 
some questions raised about the Vafa-Witten argument (see 2] for a review and references 
therein) . 

Vafa and Witten claimed that their argument should hold for any parity violating or- 
der parameter, but this might not be true in general. In fact, order parameters including 
higher order composite operators might be non-renormalizable and it is unclear how the 
argument proceeds in such . On the other hand, non-local order parameters like 

a Wilson loop around a very small compactified space direction acquire non-trivial vev's. 
Thus, one should at least restrict the type of order parameter to only local operators which 
are renormalizable, parity odd and preserve all other symmetries of the system (including 
gauge symmetry and Lorentz invariance). One of the order parameters matching all these 
requirements in 3+1 dimensions is the topological density & vap F m F ap . However, even in 
that case the argument of Vafa and Witten can be questioned {4, 5]. The aim of this note 
is to rigorously prove the Vafa-Witten theorem for this particular order parameter. 

In this case there is no doubt about the existence of a well defined free energy for any real 
value of the 9 parameter. This is guaranteed by unitarity of the theory, which is translated 
into Osterwalder-Schrader (OS) positivity of the Euclidean functional measure of the theory 

n q 

with real #-terms [6|, |7[. We recall that OS positivity is not just the standard positivity 
used in the Vafa-Witten argument. This special property excludes the possibility that 
pathological scenarios where the free energy is ill-defined [4j might occur in these theories. 
The Vafa-Witten argument can be correctly applied to show that the free energy density 
£g = Eq/V at 9 = is lower or equal than at 9 7^ 0, i.e. < £g. 

However, this fact does not mean that the theorem is already proven. One needs to 
exclude the existence of a cusp in the energy density at 9 = 0. The appearance of the 
cusp would signal the existence of a first order phase transition and spontaneous breaking 
of parity at 9 = without violation of the Vafa-Witten inequality for free energies. The 
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only extra property one needs to prove is smoothness of Sg at 9 = 0, which would make the 
existence of such a cusp not possible. In other terms, 9 = is always an absolute minimum 
of the free energy density but only if Sg is smooth we will have S' Q = 0. Only in that case we 
can make sure that the vacuum expectation value of the topological density vanishes and 
the Vafa-Witten theorem holds for this particular order parameter. 

In fact, due to Bragg symmetry S n+ g = S n -g, which follows from CP symmetry Sg = S-g 
and ^-periodicity Sg = Sg +27T , the two classically CP symmetric points 9 = 0, 7r are extremals 
of the vacuum energy density. However, Sg is not always smooth at 9 = n: there are cases 
where CP is spontaneously broken at that point. In principle, the same pathology could 
appear also at 9 = 0. 

In order to prove smoothness, the role of Lee- Yang singularities becomes very relevant. 
Because of Bragg symmetry and unitarity, non-analyticities are harder to trace in the pure 
real 9 sector. Non-analyticity in 9 can be better inferred from the lack of an analytic 
continuation into the complex 9 plane. Indeed, in that case OS positivity (i.e. unitarity) is 
not preserved and there is even no guarantee that the free energy is well defined, in which 
case it will not make sense to talk about smoothness. Essentially, there are two possible 
ways in which the 9 theory might give rise to pathological non-analyticities in the complex 
sector. One possibility is that the analytic continuation to complex values of 9 is not defined 
at all, i.e. the partition function itself becomes divergent. We will see that this possibility 
does not occur in gauge theories. The other possibility is that the partition function exhibits 
a sequence of zeros converging to 9 = in the infinite volume limit. This is the scenario 
advocated by Lee and Yang 8| which signals in many cases the existence of a phase transition 
like in the Ising model in presence of an external (imaginary) magnetic field. 

The partition function splits into a sum of functional integrals over the different topolog- 
ical sectors 

Z e (g) = J2 e- qlm9+iqRcd / 8 A e ~ s ™^ 2 3 2 det($ + m) (1) 

q= -°° c 2 (A)=q 

In a strict sense, the functional integral is UV divergent but it can be regularized in such 
a way that its positivity properties are preserved [9(. We will consider such a geometrical 
regularization whenever it is required. On the other hand, in the infinite volume limit Zg(g) 
vanishes unless the vacuum energy is renormalized to 0. But it is just the dependence on 9 



of the vacuum energy what we want to analyze. Therefore, we will consider throughout the 
paper a compact space-time with large but finite volume VT < oo. 

Positivity of the Yang-Mills Euclidean measure 0, [l(3| and fermionic determinants [1] 
implies the existence of a real effective action S e fi(g,q) = Sym(q) /2g 2 — logdet($ + m). 
On the other hand the Yang-Mills action Sym(<?) is bounded below by the BPS bound 
Syyi{q) > 87r 2 |g| in each g-topological sector. We shall see that the total effective action 
S e s(g, q) is also bounded from below by a similar linear bound involving the absolute value 
of the topological charge \q\. In the regularized theory, the coefficient of this linear bound 
increases with the UV regulating scale. This growth of the effective action is enhanced in the 
chiral limit of light quark masses by the contribution of fermionic determinants. In general, 
from positivity of fermionic determinants and Yang-Mills measure and the BPS bounds we 
have the inequality 



\Z e {g)\< elqlm(>l / SAe- s ^ {9 ' q) < / 5A e _Scff(5 ' 9) = Z (g) (2) 

9= -°° c 2 (A)=q 1=-°° c 2 {A)=q 

where 1/g 2 = 1/g 2 — (l/87r 2 )| Im#|. This shift of the gauge coupling constant can be 
considered as a change of the renormalization scale. This shows that the partition function 
of the original theory with a complex #-term is bounded by the partition function of a 
similar theory with 9= but with a different coupling constant. Since the theory at 9 = 
is unitary and renormalizable, its partition function Zq is finite and from (j2J) it follows that 
\Zq\ < oo is also finite for small values of Im#. If in addition the theory is asymptotically 
free, then the renormalization of g 2 in the UV fixed point can absorb any value of Im 6, which 
implies that the radius of convergence of the sum ([T]) in the complex 9 plane, in fact, becomes 
infinity. To make the argument more precise one should consider the regularized theory. The 
analyticity of Zq also appears naturally in lattice regularizations, even with non-topological 
regularizations of the #-term. However, to prove of absence of poles the topological charge 
bound (T5]) is essential. The bare coupling constant g has to be fine tuned according to the 
renormalization group to yield the appropriate continuum limit. But because of asymptotic 
freedom this coupling goes to zero as the UV regulator is removed and the shift from g to 
g induced by Im 9 simply implies a change in the effective scale of the continuum theory. 
However, the proof of the Vafa-Witten theorem only requires a finite radius of convergence of 
Zq around 9 = 0. This feature does not require asymptotic freedom, just renormalizability. 



We remark that one cannot use the same kind of arguments for other order parame- 
ters. For instance, the partition function associated to higher odd powers of the topological 
density might diverge for any complex value of the perturbation parameter, as they are 
non-renormalizable. A similar remark applies to the perturbation by a chemical potential 
term in theories at finite density. In that case the missing property is the lack of a BPS 
bound for the perturbation. 

The only remaining possible source of non-analyticity is the presence of Lee- Yang singu- 
larities, i.e. zeros in the partition function Zq which could prevent the existence of a unique 
limit of the free energy logZg or its derivatives at 9 = 0. 

Because of Bragg symmetry and ^-periodicity Zq is a function of sin 2 ($/2) and its zeros 
come in 4-plets {d n) — $ n , — ■#*} in the modular strip $ e (— ir, ir] x IR) of the complex 
plane. They can be parametrized by their representatives in the positive quadrant d n , 
which cannot lie on the real or imaginary axes for any finite volume. 

Only zeros converging to the origin in the thermodynamic limit can give rise to a cusp. 
Moreover, only an infinity of such zeros can produce the cusp. 

The Weierstrass factorization theorem implies that the vacuum energy density 

£v(0) = -i]nZv(0), (3) 



can be rewritten in the functional form 
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M*) = ~ | Re E P» { ln ( X ~ + Mshl2 W2) ' V) ] + gv(rin 2 (fl/2)), (4) 

where d n denotes the order n Lee- Yang zero in the positive quadrant, p n its degeneracy and 
h n and g are, respectively, polynomial and analytic functions of sin 2 (i?/2). The polynomial 
functions h n are necessary to compensate the divergent effect of zeros with slow approach 
to the infinite modulus for a given volume. 

In the thermodynamic limit the degeneracy factor p n of zeros converging to the origin 
defines a density of zeros p(z) as a function of the continuous variable z = $ n /V. Then, in 
the thermodynamic limit 

6(0) = -2Re jT J™dz 1 dz 2 p(sm 2 (z/2)) jln (l - + h(sm 2 ($/2), z)\ 

+ s(sin 2 (tf/2)). (5) 
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Notice that only if p(0) 7^ the contribution from the logarithm at z = could give 
rise to a cusp at i? = 0. The rest of the contributions from the integral have continuous 
derivatives for real 1? at 1? = 0. In fact, the logarithmic term in (jSJ) corresponds to a 

re cusp 



mm, 



Coulomb potential, in two dimensions, generated by the charge density —2p(z). T 
at 9 = is only possible if there is a non-vanishing linear charge density at 9 = 
i.e. p{z) ~ Pq8{z\ — mz§) for \z\ <C 1 with p ^ 0. A cusp corresponds to a gap in the 
potential derivative, which is proportional to the charge density — 2po- Since this charge 
density is non-positive, the corresponding potential is always repulsive, which means that 
any cusp is always a A cusp instead of a V cusp. Since the vacuum energy density £("&) 
has a local minimum at 9 = because of the Vafa-Witten inequality such a cusp cannot 
exist. The theorem follows from two very simple facts, the Vafa-Witten inequality and the 
repulsive character of the negative charge density defined by —2p(z), which follows from 
the fact that the partition function Zy($) is analytic and has zeros but no poles. The only 
way of having a a V cusp at d = is with an attractive charge density which can only 
only appear from poles in In fact, this is what happens in 1+1 massless scalar field 

theories with pseudo-periodic boundary conditions [ijj where the partition function is non- 
analytic in the complex plane. However, in QCD the analyticity of the partition function 
is based in very simple fundamental properties of the theory: unitarity, renormalizability, 
positivity and existence of BPS bounds. This makes impossible the existence of a cusp in 
the vacuum energy density at 1? = 0. In summary, the theory cannot undergo a first order 
phase transition at 9 = with parity symmetry breaking and we have a complete proof of 
the Vafa-Witten theorem. 

From these arguments we cannot, however, exclude the existence of first order phase 
transitions for other values of 9 7^ 0. In particular, we cannot discard CP symmetry breaking 
for 9 = 71. 

We remark that many other theories share the same properties and therefore the Vafa- 
Witten theorem can also be extended for all of them. One particularly interesting case is 
that of CP N non-linear sigma models. Indeed, from a similar analysis one can conclude that 
there is no spontaneous CP symmetry breaking in CP N models at 9 = 0. This result can be 
checked for the CP 1 model where the exact solution is known for 9 = [if], \v\ . 

The solution of the CP 1 model at 9 = describes the interaction of a massive scalar 
particle in the adjoint representation of SU(2), and CP symmetry is preserved confirming 
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the result of the Vafa-Witten theorem. 

Another case where an exact solution is also known is the CP N model in the large N 
limit 181 ] . In this case the system also describes a weakly interacting and parity preserving 
massive scalar particle in the adjoint representation of SU(N). Let us analyse in some detail 
how the theorem works in this particular case. The effective action in a topological sector 
q 7^ with periodic boundary conditions, is given at leading order in N by 



Nm 2 V Airq 1 [2 + IggJ Nm 2 V 
S cS (m, q) = log ^ - giVlog — -j= — (6) 



where 2m is the mass of the physical particles, /i = m exp(2ir/g 2 ) is an energy scale fixing 
the value of the coupling constant g, V = LT is the volume of the Euclidean space-time 
T 2 and T is the Euler gamma function. The effective action has a simple linear asymptotic 
regime for large values of the topological charge q > m 2 V, 

r, , \ Nm 2 V , 47r|g| N\q\ . „ ,„ N 

Seff m, q) = log + log 2 + 0(1) 7 

Air \i V 2 

The radius of convergence of the partition function Zg in the pure imaginary axis 8 = ii9 is 
i? c = A^/21og2. The critical value 6 C increases with N, and in the large N limit the whole 
pure imaginary axis is free of Lee- Yang singularities in agreement with the general theorem. 
In the complex #-plane the radius of convergence around the origin is larger than n/A even 
in the infinite volume limit, which also confirms the Vafa-Witten theorem in the large N 
limit. Now, the theorem holds not only in that limit, but for any finite value of iV as well. 

Finally, we remark that the asymptotic behavior of the effective action for small values 
of the topological charge q < m 2 V is quadratic instead of linear in q, 
„ , \ Nm 2 V , m 2 . Nnq 2 „. 4 . 

S off (m,,) = — (l-lo g -) + g -^ + Otf). (8) 

Since V = LT the change of asymptotic behavior from <£7§ to (jHJ) can be associated with 
a finite temperature /3 = l/T crossover from the low temperature regime (3 > m 2 L/\q\ to 
the high temperature regime (3 < m 2 L/\q\ and cannot be related to any phase transition 
One might expect a similar behavior in QCD, although in that case there is a finite 
temperature phase transition [20! ] . 

The arguments used above cannot exclude the existence of higher-order phase transi- 
tions which do not require spontaneous parity symmetry breaking. Indeed, if p = but 



p(z) ~ pi z% S(zi — mz 1 ^) with < a < 1 for \z\ <C 1, the vacuum energy density presents 
a second order singularity and the topological susceptibility diverges at d = 0. This be- 
havior is compatible with the Vafa-Witten inequality and might explain the behavior of the 
topological susceptibility of the CP 1 model in numerical simulations 21]. 

We thank J. Barbon, I. Cirac, V. Laliena, E. Seiler and E. Witten for illuminating dis- 
cussions. 
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